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Based on the precision experimental data of energy-level differences in hydrogenlike 
atoms, especially the IS' — 2S transition of hydrogen and deuterium, the necessity 
of introducing a reduced Dirac equation with reduced mass as the substitution of 
original electron mass is stressed. Based on new cognition about the essence of 
special relativity, we provide a reasonable argument for reduced Dirac equation to 
have two symmetries, the invariance under the (newly defined) space-time inversion 
and that under the pure space inversion, in a noninertial frame. By using reduced 
Dirac equation and within the framework of quantum electrodynamics in covariant 
form, the Lamb shift can be evaluated (at one-loop level) as the radiative correction 
on a bound electron staying in an off-mass-shell state-a new approach eliminating 
the infrared divergence. Hence the whole calculation, though with limited accuracy, 
is simplified, getting rid of all divergences and free of ambiguity. 
Keywords: Reduced Dirac Equation, Lamb shift, off-mass-shell 
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I. INTRODUCTION 

As is well known, the Dirac equation for electron in a hydrogenlike atom is usually treated 
as a one-body equation with the nucleus being an inert core having infinite mass and exerting 
a potential V{r) = — ^ {h = c = 1) on the electron. Then the rigorous solution of energy 
levels reads [1]: 
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where j is the total angular momentum. The expansion of f{n,j) to the power of {ZaY 
is given as P 
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Obviously, besides the rest energy of the electron given by the first term, the second term 
has exactly the form of Bohr energy level except that the mass rrie must be replaced by the 
reduced mass 
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rUe + rriN M 
with rriN being the mass of the nucleus and M = rrie + mjv. 

However, as discussed in Refs.pLj and |2j, the concept of reduced mass in relativistic 
quantum mechanics (RQM) is ambiguous to some extent. Beginning from 1950's, a number 
of authors have been devoting a great effort at the level of two-body RQM and that of 
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quantum electrodynamics (QED) to take account of the recoil effect [HOIl], incorporating 
their results in a compact form (to order of a^): 



E = M + f,[f{n,j) -1] - ^\f{n,j) - + 
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A comprehensive review on the theory of hydrogenlike atoms can be found in Ref.[2Zl- The 
aim of this paper is two-fold: First, based on the experimental data of hydrogen IS — 2S 
transition frequency [5] and its isotope shift of hydrogen and deuterium [6], we stress the 
necessity of the introduction of reduced mass /i (section II) before we are able to argue 
the reasonableness of introducing a "reduced Dirac equation" with /i as the substitution of 
rrie (section III). Second, based on above conception, we will present a calculation of Lamb 
Shift (LS) as an off-mass-shell effect by performing the evaluation of self-energy diagrams of 
electron (section IV) and photon (section V) as well as the vertex function (section VI) at 
the one-loop level of QED in covariant form. The new insight of our calculation is focused 
on the regularization renormalization method (RRM). As initiated by J-F Yang j[7j and 
elaborated in a series of papers (P [2H - [26| ISaf l8b] and references therein), we can get rid of 
all ultra violet divergences in the calculation of quantum field theory (QFT). Furthermore, 
in this paper, we will be able to get rid of the annoying infrared divergence in the vertex 
function by treating the electron moving off its mass-shell to certain extent which is fixed 
through the evaluation of self-energy diagram or by the Virial theorem. Based on above 
improvements, the one-loop calculation yields values of LS in a simple but semi-quantitative 
way (section VII and VIII). Although the accuracy is limited at one-loop level, we hope our 
approach could be served as a new starting point for calculations at high-loop order to get 
accurate results at a comparably low labor cost. The final section IX and Appendix will 
contain a summary and discussion. 



II. THE IS" - 25 TRANSITION OF ATOMIC HYDROGEN AND DEUTERIUM 

In the last decade, thanks to remarkable advances in high resolution laser spectroscopy 
and optical frequency metrology, the IS* — 2S two-photon transition in atomic hydrogen H 
(or deuterium D) with its natural linewidth of only 1.3Hz had been measured to a very 
high precision. In 1997, Udem et a/.determined the IS' — 25" interval of H being [5] 
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f^\lS - 2S) = 2466061413187.34(84) kHz (7) 

Even four years earlier, Schmidt-Kalar et al. measured the isotope-shift of the IS — 2S 
transition of H and D to an accuracy of 3.7 x 10~*[B], giving (as quoted in flU]): 

A/ = /(^)(25- 15) -/(^)(25- 15) = 670994337(22) kHz (8) 

(In 1998, Huber et a/.measured a more accurate data [28]: 670994334.64(15) kHz), which 
is of the order of 10~^ in comparison with Eq.([7|. As pointed out in Ref. [6j, this 671 GHz 
isotope-shift can be ascribed almost entirely to the different masses of proton (p) and 
deuteron (rf). And the nuclear volume effects become important because the QED effects 
cancel considerably in the isotope shift. 

Here, we wish to emphasize that in the first approximation, both experimental data ([T]) 
and ([s]) can be well accounted for by simply resorting to Eq.([T| with rUe replaced by the 
reduced mass 



l^H = ^ , f^D = ^ (9) 

rrie + Trip me + md 



for H and D respectively. 

Indeed, adopting the following updated values p^HT3] 
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0.51099906 MeV 
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= 0.532513542 x 10"* 
a'^ = 0.151005223 x 10" 
: 1.2355897 x 10^° Hz 
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3.28984124 x 10'^ Hz 



-- 1836.1526665 

bn = 5.446170255 x 10-^ 
bo = 2.724436319 x 10"^ 
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0.999455679 
0.99972763 
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we can calculate the energy difference of 25" and IS* of H through Eq.([T]) with rUe replaced 
by fJ^H (the superscript RDE refers to the reduced Dirac equation) 
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Ai?^^^(25 - 15) = /i^,[/(2,l/2)-/(l,l/2)] 

= ^(1.996950159 x 10"^) 

= 1.2355897 x 10^° x 0.999455679 x 1.996950159 x 10"^ 

= 2.466067984 x 10^^ Hz (17) 

which is only a bit larger than the experimental data Eq.([7]) with accuracy 3 x 10~^. 
However, a more stringent test of RDE should be the isotope shift of H and D. We have 



(bn - bo) - {bjj - bjj) + {bjj - bj^) + ■ ■ ■ = 2.719511528 x 10"^ (1 



l + bo l + bn 



A^^^f = {fxn - /iH)[/(2, 1/2) - /(1, 1/2)] = 6.7101527879 x lO^^ Hz (19) 



which has only a discrepancy larger than the experimental data, Eq.(|8j) by 20.941 MHz 
with accuracy 3 x 10~^. Of course, it is still not satisfied in an analysis of high precision [6]. 
Let us resort to the Eq.(6), where the third term does provide a further modification: 

= ^me[(6H - bo) - 3{bjj - bl) + ■ ■ ■ ](-6.646361554 x 10"^°) = -11.176 MHz (20) 

which brings the discrepancy between the theory and experimental down to less than 
10 MHz. 

Although the detail explanation for this discrepancy remains quite complicated [6j, the 
above comparison is enough to convince us that the inevitable appearance of reduced mass 
in the RDE or Eq.(6) is by no means a simple fortune. It must have a deep reason from a 
theoretical point of view. Notice further that once the conditions rrie ^ mp and rrie ^ rrid 
hold, the difference of spin between p and d seems not so important. So in next section, we 
will strive to justify the reduced Dirac equation on a reasonable basis. Of course, it is still 
an approximate one, but seems much better than the original Dirac equation when dealing 
with hydrogenlike atoms. 
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III. REDUCED MASS AND REDUCED DIRAC EQUATION 

Consider a system of two particles with rest masses mi and Their coordinates in 
the center-of-mass (CM) system are ri and r2 respectively, as shown in Fig.l. If there 
is a potential V{r) = V{\ri — r2|) between them, two equations miri = — Vr^(r) and 
m2r2 = VrV{r) will reduce to one: 

A'|^ = -V.m (21) 

at'' nil + ^2 

At first sight, the definition of center-of -mass (CM) in classical mechanics rriiri = m2r2 
becomes doubtful in the theory of special relativity (SR) because the mass is no longer 
a constant. But actually, we can still introduce the coordinate of CM in the laboratory 
coordinate system (LCS) (with r'^ and r2 being the coordinates of mi and m2): 

R= ^{mir[ + m2r'^) = {X,Y,Z), (M = mi + m2) (22) 
and the relative coordinate of mi and m2 (r^ = — R, i = 1, 2) : 

r = r'l - r's = ri - r2 = (x, y, z) (23) 
Here the motion of CM in the LCS is assumed to be slow and so 

d mi Odd m2 d d , , 



dx'i M dX dx' dx'^ M dX dx 
Notice that the momentum P of CM and the relative momentum becomes operator 

in quantum mechanics (QM) without explicit dependence on mass: 

p = -iWa, p,. = -iWr (25) 
Thus the momenta of mi and m2 in laboratory coordinate system (LCS) read: 

p; = -iWr; = ^P + Pr, P2 = -iWr^ = ^P - p,. (26) 

Since the center-of-mass coordinate system (CMCS) is also an inertial frame which can be 
transformed from the LCS via a linear Lorentz transformation, it is defined by the condition 
that P = in CMCS. In other words, CMCS is defined by the condition pi + P2 = 0, or 



from Eq. ( 26 ) : 
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Pi = -ihVr[ = Pr, P2 = -«Wr!, = -Pr (27) 

Evidently, the above definition of CMCS remains valid in the realm of relativistic QM 
(RQM) even the exact meaning of CM seems obscure to some extent due to the conjugation 
relation of a particle's position and its momentum, see Fig.l. 



Now, from Eq.(27), it is natural to replace pi and p2 by p^, reducing the two-particle 



degrees of freedom to one. In the meantime, the origin of CMCS is discarded, it is substituted 
by the position of m2 (r = ri — r2). We will call the system associated with r the relative 
motion coordinate system (RMCS), which should be viewed as a deformation of CMCS. The 
transformation from CMCS to RMCS is by no means a linear one. Rather, the origin of 
RMCS {1712) is moving non- uniformly in the CMCS. Therefore, while rest masses mi and 
1712 remain the same in both LCS and CMCS, they reduce to one mass fi = ^'^^ for mi 
in RMCS (or for 1712 if mi is chosen as the origin of RMCS). 

Let us express the total energy E = Ei + E2 in CMCS in terms of Pr and reduced mass 



(/i = M = mi+ ms), where 
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(28) 



Treating all pi,P2 and Pr being c- numbers, we have 



E^ = {El + E2f = M^ + —pI + —,pt{A ) + 



(29) 



where the expansion in p^ is kept to the order of p^. Two extreme cases will be considered 
separately: 

A. m2 > mi, p < mi, M > /x: 



E^ 
E 



M 



14 

1 + 
1 



1 



1 o 1 



p4Q _ ^) + 



rPr - 



Pta 



M' 



+ ■■■ 



8/i' 



M - 11+ ^/n'^ +Pr ^m2 + {mi - n) + -x//^^ + Pr 



(30) 



E' = E - m2 = {mi - n) + \fW^-^ 



(31) 



8 



B. mi = m2 = m, = y, M = 2m = 4yU Then to the accuracy of p^, we have : 



= + —pi = Am^ + Apl 



,2 ^ „4 ^ 1 „2 



£'.£_M=-rf-— (ifrf«,^) (32) 

It is interesting to see that after introducing and pr, the energy E' in RMCS looks quite 
"relativistic" in the case A whereas it looks rather "non-relativistic" in the case B even both 



of them are derived from the relativistic expressions, Eq.(28), approximately. 



Since the RMCS is not an inertial system, the original mass of mi in CM changes abruptly 



to /i as shown in Eq.(31). How can we derive the reduced Dirac equation (RDE) in RMCS? 
Fortunately, we already found a basic symmetry, the space-time inversion symmetry, which 
not only serves as the essence of special relativity (SR), but also goes beyond it to derive 
the original Dirac equation and the tachyon theory for neutrinos [HHTTj. Based on this 
symmetry, we are going to derive the equation in RQM for case either A or B respectively. 

Let us consider case B (mi ~ 1712) first. The motivation is stemming from the success 
of using the Schrodinger equation to heavy-quarkoniums like cc and hh in particle physics 
(|18j. see also [15] §9.5 D). Ignoring the spin of both mi and m2, we assume the coupling 
equations in laboratory system for the two-particle system as: 



^h% = (mi + m2)cV + V{K " r'Mv + x) - + StXO^^ + x) ^33. 

= -(mi + m,yx - V{K - r',\){v + x) + {£-X[ + + ^) 

where = (p{r[,r2,t) and x = x(i'i)i'2)^) hidden "particle" and " antiparticle" fields 
of the two-particle system (From now on, the r^(z = 1, 2) is the flowing coordinate of "fields" 



in QM, i.e., that of "fictitious point particles". See Fig.l). Eq.(33) remains invariant under 



the (newly defined) space-time inversion (r'^ — t- —r[,r2 — > — r'gjt — t- —t): 

y?(-r;,-r'2,-t) — > x{r[,r'2,t) ^^^^ 
x{-r\, -r'2, -t) — y ip{r[,r'2,t) 

V{-r[,-v'2,-t)^V{v[,r'2,t) (35) 
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Note that, however, the time t is not contained in V exphcitly. Eq.(35) merely means 



that both nil and 1712 {mi ^ 1712) transform into their antiparticles under the space-time 
inversion. Actually, the hidden antiparticle field x enhances in nearly equal strength in 
nil and m2 when fictitious particles' velocities increase with the enhancement of attractive 
potential V{r). 

After introducing the CM coordinate R = -^(mir'i +m2r'2), (M = mi +7712) and relative 
coordinate r = r'^^ — r'2, and setting 



93 = $ + i 



n 



X = '^ 



t 



(36) 



we find (/i - 



M 



jji ft 



(37) 



Its stationary solution reads 



$(R, r,t) = V'(r)exp 



h 



;p ■ R - Et) 



(38) 



where E is the total energy of the system while P the momentum of CM. The reduced 
"one-body" equation for ip{r) turns out to be: [*] 



'2/i 



Vl + V{r) 



(^2 _ ^2^4 _ p2^2^ 



(39) 



2Mc2 



We set P = (z.e.turn to CMCS) and denote the binding energy B = Mc^ — E, yielding: 



B = Mc^ 



I- a 



2e 



a/2 



-e + 



1 



(40) 



[*] With Eq. ( 35 ) , Eq. ( [37| ) is invariant under the space-time inversion (r 
the mass inversion (mi — > — mi,TO2 ^ 



-r,t — >■ —t). Equivalently, under 



7712), Eq.pTt and Eq.(39l remain invariant in the sense that 



not only /i — > —^,M — > —A/, but also Vir) — > — V^(r),e — > — e. Notice that, however, the simultaneous 
inversion of mi and m2 implies mi ~ m2, so both particles change under their mutual interaction V{r) 
simultaneously. Here V , being the "internal potential energy" of two-body system, was called as a "scalar 
potential" . We see that either the invariance under the space-time inversion or that under the mass 
inversion is capable of showing the particlc-antiparticle symmetry (i.e., relativistic nature) of a system 
essentially. 
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Notice that although Eq.(39) looks like a "non-relativistic" stationary Schrodinger equa- 



tion, it is essentially relativistic. This can be seen from its remarkable property that the 
eigenvalue e has a lower bound — ^Mc^, corresponding to E^^^ = (-Bmax = M)\ 

An example is: consider "positronium" composed of e"*" and e~ with charge Ze and 
—Ze respectively. Once when the "fictitious charge number" Z increases from 1 to Z^ax = 
{-^y^'^ = 16.555, the whole bound system would have lowest ground energy Emin = 0! So 



Eq.(39) is really a relativistic QM equation capable of giving a nonperturbative solution 



under the strong coupling. 



Eq.(39) provides a justification (realization) of conjecture Eq.(32) relevant to case B 



(mi ^ 1712) where the spin of both particles is merely of second importance. 

Now let us turn to case A where 7712 ^ mi, taking the spin of mi into account but 
ignoring that of 1712 as before. Based on the experience in case B, also because of great 
difficulty to derive the equation starting from the laboratory system for this case A, we 
directly introduce the reduced Dirac equation (RDE) in the RMCS as a pair of coupled 
equations of two-component spinors (p{r,t) and x(r, t), {c = h = 1) 



(41) 



iif = i(7i ■ VrX + /^V^ + ^(r)v' 
iX = iai ■ VrV? - I^X + V{r)x 
with /i replacing mi. Here ui are Pauli matrices acting on the spin space of particle 



mi. Eq.(41) is invariant under the space-time inversion (r — t- —r,t — t- — t), ip{—r,—t) — )■ 
x(r,t), x(— r, — t) —7- '^{r,t) whereas we assume 



y(-r,-t) 



-V{r,t) 



(42) 



here in contrast to Eq.(35) for the case B. [f] 



The reasons are as follows: (a) Eq.(41 ) should degenerate into the original Dirac equation 
when 1712 00, fj, = > nil. (b) Since now m2 ^ rrii (but 1712 7^ 00), mi is moving 



For a hydrogenlike atom, V{r) — — does not contain time t explicitly. Eq. ( 42 ) merely means that under 
the space-time inversion, the electron transforms into a position whereas the nucleus remains unchanged. 



See point (a) of section IX. Previously, the V in Eq.(42) was called as a "vector potential", meaning the 



"potential energy" of the electron in an "external field" of nucleus. Note that, formally, Eq.(41 1 remains 
invariant under a mass inversion as /i - 
frame RMCS. Actually, since rrii — m 



fi,(j) XjX ~^ 4' (^(r) remains unchanged) in a noninertial 
—rrie, but 1712 



2me ■ 
M ' 



So Eq.(41) 



has an inaccuracy up to (< 1.1 x 10 ^ for H) 
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much faster than m2 in the CMCS. Hence the antiparticle field x enhances much appreciably 
in mi than that in m2, a situation totally different from that in the case B where mi ^ m2. 



c) If instead of Eq.(42), we still assume V^(— r, —t) — )■ V(r,t) like Eq.(35) and change the 



sign before V^(r) in the second equation of Eq.(41) to keep its invariance under the space- 
time inversion, then we would get an equation which would lead to a reversed fine-structure 
of atom {e.g., the P1/2 state would lie above the P3/2 state), a wrong prediction obviously 
excluded by experiments. 

However, one kind of invariance is not enough to fix an equation. Indeed, the beauty of 
Dirac equation or RDE is hidden in two symmetries: besides the symmetry of space-time 
inversion, it has another left-right (parity) symmetry. To see it, we define 



V2 



and recast Eq.(41) into: 



(43) 



= icTi ■ Vr^ + /ir] + 1/(r)^ 
ii] = —iai ■ VrV + /U^ + V{r)r] 
which is invariant under a pure space inversion (r — r, t — )■ t) if assuming 



(44) 



e(-r,t)^r7(r,t), r/(-r, t) ^ ^(r, t), V (-r) ^ V (r) = V (r) (45) 

The parity invariance of Dirac equation or RDE has a far-reaching consequence that 
the Dirac particle is always a subluminal one. By contrast, once the parity is violated to 
maximum, a superluminal particle (tachyon) will emerge. Interestingly enough, any theory 
capable of treating particle and antiparticle on an equal footing must respect to the common 
basic symmetry — the invariance of space-time inversion. The new insight of this section is 
this symmetry can be applied even in a noninertial frame — the RMCS. Of course, the validity 
of RDE can only be verified by experiments as discussed in section II, although it is still an 
approximate description of nature like any other theory in physics. For further discussion, 
see section IX. 
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IV. SELF-ENERGY CORRECTION OF A BOUND ELECTRON IN ATOM 

In our understanding, one important reason why the calculations of QED for electron in a 
hydrogenlike atom is so complicated lies in the fact that while calculations are performed in 
the CMCS, the center of potential (the nucleus with mass m2 = itln) undergoes a complex 
motion. So the recoil effect interwinds with the high-loop correction of QED, as discussed 
in many chapters of the books [Ij and [2j. We will try to find an alternative approach by 
adopting the RDE and doing calculation in the RMCS. Let us begin with the Feynman 
diagram integral (FDI) of electron self-energy at one-loop level, adopting the Bjorken-Drell 
metric and rationalized Gaussian units with electron charge — e(e > 0), see Fig.2(a) ([5a]). 

- ^S(p) = {tef [ -^,^-^1^—4. (46) 

^ ' J (27r)4 ' ^- ^- /i ' ^ ' 

Here a free electron with reduced mass /i is moving at a four-dimensional momentum p, 

whose spatial component is just the relative momentum discussed in the previous section. 



k is the momentum of virtual photon. As usual, a Feynman parameter x will bring Eq.(46) 
into 

1 dx 

D=l |F-2p.t..+ (p^-,^).p - 'V = -2b'-W+4, (48) 

(/ = J9'^7^, p ■ k = p^k^). A shift in momentum integration: /c K = k — xp recast 



Eq.(47) into 



- iS(p) = -e^ f dx[-2{l - x) j/+ 4/i]/ (49) 
Jo 

with a logarithmically divergent integral (in Minkowski momentum space): 

'-jwfW^- M^^pV + (,^-/). (50) 
Our new regularization-renormalization method (RRM) is based on a cognition that the 
virtual process in the self-energy diagram does provide a radiative correction to the electron 
mass but only when the electron is off the mass shell, i.e., 7^ /i^. When it is on the mass 
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shell, p = /i , the appearance of a divergent integral like / in Eq.(50) is essentially a warning 
on the fact that to calculate the mass of electron is beyond the ability of perturbative QED. 

Let us consider the converse: if does modify the electron mass /i to some extent, it 
must comes from the divergent integral /. However, the latter is a dimensionless number, 
we can change the unit of M (and k) at our disposal without any change in the value of 
/. So any real change of fi (on the mass shell) is incredible. The deeper reason lies in a 
"principle of relativity" in epistemology: everything is moving and becomes recognizable 
only in relationship with other things. What we can understand is either no mass scale or 
two mass scales, but never one mass scale. For instance, in the famous Gross-Neveu model 
[19] . a massive fermion is created only in accompanying with the change (phase transition) of 
its environment (vacuum) which provides another mass scale (a standard weight). Another 
example is just the change of electron mass from rrie to fin in a hydrogen atom due to the 
coexistence of atom nucleus — the proton, this change is also a nonperturbative effect. 



Therefore, we expected too much in the past. There is no way to evaluate Eq.(50) 



unambiguously or pick out some finite and fixed modification on the mass /i. What we can 



do is to separate the valuable information carried by Eq.(50) from an arbitrary constant 
which will be introduced by a simple trick and then fixed by the experimental data of fi. 
We will see the information telling us exactly how the value of I changes when the electron 
is moving off the mass shell. 



To handle Eq.(50 ), we perform a differentiation with respect to the mass-square parameter 
M^, then the integration with respect to K becomes convergent, yielding: 

— = — — (51) 

which tells us that while the exact value of I remains obscure, its change linked with 
has a definite meaning. So we reintegrate Eq.(51) with respect to and arrive at 



—7 —7 

/ = 7^(lnM2 + Ci) = In — (52) 

(47r)2^ (4vr)2 fxj 

where an arbitrary constant Ci = — In fi^ is introduced (/i2 should not be confused with 
the reduced mass n). Further integration with respect to Feynman parameter x leads to 



14 



S(p) 


= A + Bj/ 










A 




2 - 21n 






In 




B 


a 

47r 


21n^ - 


3 


p2 


1 + (^^r ) In ^'^'-f^ 



(53) 



Using the chain approximation, we can derive the modification of electron propagator as 

i i 1 iZ2 



where 



1 



(54) 



' 1-5 

is the renormahzation factor for wave function of electron and 



(55) 



11 + A 



1 - B 

is the renormalized mass of ii. The increment of mass reads 



(56) 



6n = fiR- ij, 



A + jiB 
1 - B 



(57) 



For a free electron (in the atom), the mass-shell condition = /x^ should lead to 

5uL2_„2 = —(5 — 6 In — ) = 

as discussed above [J] . So we must set /i2 = fj.e~^^^ which in turn fixes 

1 



■^2 L2 = u2 



1 + 



1 



a 

3^ 



(58) 



(59) 



However, the above evaluation further provides us with important knowledge of 6fi when 
electron is moving off the mass-shell. Consider the similar diagram in Fig. (2b), we can set 
on an average meaning that 



(60) 



[I] We will keep the same mass symbol /i through out high-loop calculations of QED and reconfirm (renor- 
malize) it at every step by experiment. Just like one has to reconfirm his plane ticket before his departure 
from the airport, he must use the same name through out his entire journey |8bj . 



15 



with C > 0, which imphes from Eq.(57) with Eq.(53) that [20] : 



(61) 



An 1 + a/3n 

where some terms of the order of ('^ or C^l^C ^ire neglected since C ^ 1. Eq.(61) estab 



hshes the correspondence between the mass modification 6fi and the parameter ( describing 
the off-mass-shell extent of electron in the bound state. For a hydrogenlike atom, we may 
ascribe 5/i to the (minus) binding energy of electron in the Bohr theory: 



2n2 



6fi 



(62) 



Then Eq.(61) gives the value of C for fixed values Z and n. We will see from the vertex 



function that these values of ( are crucial to the calculation of Lamb shift (sections VII and 
VIII). 

V. PHOTON SELF-ENERGY 

As discussed in various text books [2T1423] , we encounter the FDI of vacuum polarization 
Fig. 2(c) as [8a]: 

n,.(g) = -(-ze)2Tr / i%i,^l.r^ (63) 

Here q is the momentum transfer along the photon line and m the mass of electron. 
Introducing the Feynman parameter x as in previous section and performing a shift in 
momentum integration: p — )■ i^T = p + xg, we get: 



n^^(g) = -4e2 I dxih + h 
Jo 



(64) 



where 



d^K 2K^K^ 



(27r)4 (7^2 _ ]Vf2)2 



(65) 



with 



= + g^(x^ - x) 



(66) 
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is quadratically divergent while 



{K^ - M2)2 



(67) 



is only logarithmically divergent like that in Eq.(50). An elegant way to handle Ji, 
Eq.(65), is modifying to 



(68) 



and differentiating Ji with respect to a. After integration with respect to we reinte- 
grate it with respect to a twice, arriving at the limit o" — )■ 0: 



h 



{Any 



[m + g (x — x)] In ■ 



^1 



(69) 



with two arbitrary constant: Ci = — ln/i| and C2. Combining Ii and I2 together, we 
find: 



T-r / \ 8^'^ / 2\ f 1 r 2 m? + q^ix^ — x) Aie^ 

^t^Aq) = lY^W^ - 9f^uq ) / dx{x - X) In 2 7^^^Mi.<-2 (70) 

The continuity equation of current induced in the vacuum polarization 



q'-U.M = (71) 

is ensured by the factor {q^iqu — g^ivq^)- So we set C2 = 0. Consider the scattering between 
two electrons via the exchange of a photon with momentum transfer g — )■ 0. Adding the 
contribution of 11^,^ (g) to the tree diagram amounts to modify the charge square: 

e'-^el = Z,e\ ^3 = 1 + f (In ^ - + ■ ■ • ) (72) 

As in Ref. [8b], we will set /is = m so that at the Thomson limit :limq_!.o e\ = e^. However, 
for the purpose of calculating Lamb shift (LS) below, the second term in the parenthesis of 
Z3 is important because for a bound state it contributes a term of effective potential (adding 
to Coulomb potential), called the Uehling potential ([23j.p.253): 

5(r) (73) 



15m 
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VI. THE OFF-MASS-SHELL VERTEX FUNCTION 

Consider an electron (see Fig. 2(d)) moving in a hydrogen atom, its momentum changes 
from p to p' via the scattering by the proton and an exchange of virtual photon with mo- 
mentum k. The FDI at one-loop level reads 

A,ip,p) = Her J (74) 



However, different from |8b] and many other literatures, not only the reduced mass fi 
(instead of m) of electron is used, but also a new approach will be adopted. We assume that 
the electron is moving off- mass-shell in the sense of (as in section IV): 

p'=p" = Al~0 (75) 

We still have 

P'-P = (l, P-(l = -\^(f (76) 

Introducing Feynman parameters u = x + y and v = x — y, we perform a shift in the 
momentum integration: A; K = k — {p + q/2)u — {q/2)v, thus 

A^ = -te^[h^^ + h] (77) 



= [Al - - ^j]u' + ^jv' + Cfi'u (79) 
h = du dv (80) 



(4-4m-2m2);,2 ^ 1-0^^ + 22 (M2-M)/ig'^a,,, - (2 - 2m + — - -)g2^^ - (2 +2M)t;/ig^ (81) 



Set K"^ = K"^ — + M^, then I3 = 1'^ — 7^ and J3 is only logarithmically divergent 
and so can be treated as in previous sections, yielding: 
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with /ii an arbitrary constant. 

However, unlike Ref . ISbl where the calculation was conducted on the mass-shell, now the 



off- mass-shell integration in Eq.(82) can be performed in the approximation that ^ <^ 1 



and C ^ 1 (Q^ = ^9^) Q is the three-dimensional momentum transfer) which will be enough 
to calculate the Lamb shift (LS). Denoting 

a=[^^\l-0 + ^]u' + C^^'u, b=^ (83) 
we will perform the integration with respect to v and u rigorously: 

fJn 1 + ,/n2 /Any 

/ dvln(a-bv^) = 2u[lnfi^ + lnu + ln[(l-C)u + C]-4u + 2.hZln " (84) 

J-u \ Q 1 - ^/Q^/Aau 

Expanding the last term and keeping only up to the order of ( and Q'^/4fi^, we obtain 

[ du [ dv\n{a-bv^) -Infi^ -1 + C+ ^{1-0 (85) 
Jo J-u 6/i 

To our great pleasure, throughout the evaluation of J4, there is no any infrared divergence 
which would appear in previous literatures when integrating with respect to u with lower 
limit zero. To avoid the infrared divergence, e.g., in |8bj . a cutoff was introduced at the 
lower limit. Now the infrared divergence disappears due to the existence of off-mass-shell 
parameter (. For example, we encounter the following integral, in which no cutoff is needed 

(A = (l-C) + gW~l): 

^ In ^ (86) 



lo u + C/X X X 
Hence, after elementary but tedious calculation, we find: 



Mp', p) = f - In ^ - 3C + 4(1 + C) In C]7. + f ?.7.(i + K + 1 1^ C + 2C In C) ^^^^ 



+z^a-^^^(l + 3^ + 2ClnC) 



VII. CALCULATION OF LAMB SHIFT AS AN OFF-MASS-SHELL EFFECT AT 

ONE-LOOP LEVEL 



There are three parts in Eq.(87). The first part in combination with the vertex 7^ at tree 
level provides a renormalization factor as 
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= 1 



a 11 



ln^-3C + 4(l + C)lnC] 
1^1 



(88) 



Further combination with Z2 in Eq.(55) and Z3 in Eq.(72) leads to a renormahzed charge 
(at one-loop level, see Fig. 2): 



(89) 



However the Ward identity implies that [2T1423 



(90) 



Therefore 



aR = -r = ZzOi 
Ait 



(91) 



Note that Ward identity holds not only for an electron on the mass-shell, but also for 
off-mass-shell case. Hence for every bound state in hydrogenlike atom with a definite value 



of ( {Zi and Z2 are functions of (), the arbitrary constant /ii in Eq.(88) plays a flexible 



role to guarantee the validity of Eq.(90) (other two constants /i2 and /is had been fixed in 



Eq.(58) and (72) respectively). For further discussion, see section IX. 

The second part of Eq.(87) contains (5^7^. Just like the Uehling potential in Eq.(72) 



(with = —Q^), it contributes an effective potential of 6 function type as 



1 1 4 
-,[--- -C -- In C-2C In Q6{r) 



(92) 



Finally, the third part of Eq.(87) amounts to a modification of electron magnetic moment 



in the atom, the gyromagnetic ratio of electron reads: 



^7 = 2[1 



a 
2^ 



;i + 3C + 2ClnC)] 



(93) 



We will call the anomalous part of magnetic moment a = a = a{l + 3( + 2(ln(). The 
radiative correction on the magnetic moment of an electron has two consequences. One is a 
modification to the L-S coupling in a hydrogenlike atom (with charge number Z) [2T1 122] : 



H 



rad 
LS 



, a , aZ 
2{—)—^a-L 



'211' A 



(94) 
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Here the electron mass has been modified from m (see, e.g., [15]) to fi which can be 
derived from the reduced Dirac equation. 

Another consequence of anomalous magnetic moment of electron exhibits itself as an 



additional potential of 6 function type like Eq. ( 73 ) [211 122] 



Zaa 



Sir) 



(95) 



Note that Eqs.(94) and (95) are only effective to states with L 7^ and S state with 
L = respectively. 



Adding the results of Eqs.(94), (95) and the sum of Eqs.(73) and (92) multiplied by Z 



together to get all radiative corrections (at one-loop level) on electron in the hydrogenlike 
atom, then we get the effective potential as 



eff 



Zo? 



Zo? 



[-|lnC-K-2ClnC- 
;i + 3C + 2ClnC)(T-L 



1 _ 

6 15 m? 



+ i(l + 3C + 2ClnC)]5(r) 



(96) 



^.-[-|lnC + 3^ + C-ClnC]5(r) + ^(l + 3C + 2ClnC)a-L 
where we take ^ jvn? ^ 1 in the Uehling potential to make the formula simpler for a 



semi-quantitative calculation. Eq.(96) leads to the energy modification of a bound state 



(with quantum numbers n,l,j) in a hydrogenlike atom: 



6{r) ^ 1^^,(0) p 



3^,3 



(/ = 0) 



(97) 



(98) 



(99) 



_ Z^a' ^ 1 + C(3 + 2 In C) J I, (j = ^ + 1/2) 



vrn' 



3 



/(2/ + l)(/ + l) 



where 



(/ + 1), (j = /-l/2) 



(100) 



Ry = l-a^fi = —R, 
2 m 



(101) 
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VIII. ENERGY-LEVEL DIFFERENCE IN HYDROGENLIKE ATOM: THEORY 

VS. EXPERIMENT 

We will study some energy-level differences near the ground state of hydrogenlike atoms, 
where precise experimental data are available. Theoretically, the energy level is fixed pri- 
marily by the formula derived from the reduced Dirac equation (RDE), i.e., Eq.Q with mg 
substituted by /i^ where the subscript A refers to atom H, D or ife+, et al..\ 



= ^^A[f{n,J) - 1] = Y^U{n,j) - 1] 

= ^-(1.2355897 xl02°)[-^-i^(^^- ^7^(102) 
l + bA 2n2 3n3 ^ j + 1/2 An' J ^ ^ 

Further recoil corrections Eq.(|6]) derived by previous authors will be divided into two 
terms: 



2Ma 



1 



rriebA 



2fl + 6.)3 



2n^m 



N 



1 + 



t)(i-M 



(103) 
(104) 



Next comes the radiative correction calculated by QED at one-loop level, Eq.(98): 



where 



_J_g(£!^,)|(_^.„C.|.2C(i-.nC)).,„. ^^^g;^^'"^' c,(i-M 

(105) 



1 

_ 1 



J 

3 = 1 

(A) 



Finally, the finite nucleus size (NS) with radius r]^ brings a correction ^TU] : 



(106) 



(A) 



3 ^ rUe 
I 1 



3^(4.386454987 xlO^)[J^]^5,o 



(107) 



As explained in Eq.(61) with Eq.(62), the value of off-mass-shell parameter ^ in Eq.(105) 



can be calculated from the electron self-energy at one-loop level: 
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Z^a _ 1 (C<*> -2C<^>lnC 
In 1 + a/Svr 



(108) 



where the superscript < S > refers to "self-energy". However, we may derive the value 
of C in an alternative way. Divide the square average of four-dimensional momentum p into 
two parts: 



< >=< > - < > 



(109) 



where 



< >= = {^- Bf ~ /i^ - 2/x5, 



(110) 



since the binding energy 



B 



^^^^ 



(111) 



The square average of three-dimensional momentum p, < p^ >, can be evaluated by the 
Virial theorem {e.g., [U]). In a Coulomb field, an electron has potential energy V = 
and kinetic energy T = j^p'^- Then 



< p^ > = 2/i < T >= 2fi[-B- <V >]= 2iiB 



AB 
1^ 



Comparing Eq.(112) with < >= /i^(l — C"^^^)? we find 



c 



<v> 



AB 



IZ'^o? 



(112) 



(113) 



where the superscript < V > refers to "Virial theorem". Table 1 gives the values of 
C"^^^ and C^^'* with their logarithm values as well as two kinds of "average", i^<^+^> = 

^^<5>^<y>^ to be used in Eq.(|l05|). 
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Table 1. Off-mass-shell parameter and ln( 





C<^> X 10^ 


-lnC<^> 


C<^> X 106 


-lnC<^> 


^<s+y> X io5 








1 

16 


1.546093458 


8.77461 


^ = 6.6564192 


11.91992886 


8.0632 


9.425609 


3.2080284 


10.34727 


1 
4 


7.446539697 


7.20259 


= 26.6256771 


10.5336345 


38.5639 


7.860609 


14.0808 


8.86816225 


1 


37.73719345 


5.57969 


2a2 = 106.502 


9.147340142 


194.011 


6.2450103 


63.39626 


7.36351521 



Now we are in a position to discuss a number of cases: 

(a) The so-called classic Lamb shift of hydrogen atom was measured experimentally as 



L^P(25 - 2P) = Eh{2Sii2) - ^h(2Pi/2) = 1057.845 MHz (114) 
Theoretically, in this case {hn = 5.446170255 x 10-^ r% = = 0.862/m), Eqs.(|l02( 



and (103) make no contributions while Eqs.(104) and (107) only contribute 

\2Si/2 - 2P1/2) = -E^"™^'-2(2, 1/2, 1) = -2.16156 kHz 



(115) 



and 



AE^^{2S - 2P) = 0.14525347 MHz 



(116) 



respectively. The dominant contribution comes from Eq.(105). If using C^'^^j we obtain 



^^Had<5>^25 - 2P) = ^1(4.06931316 x 108)[-| lnC<^> + ^ + 3C<^> - |C<^> \nC<^>^ 



l+bH 8 



1000.6567 MHz 



If we use another three values of In ( in Table 1 , we get 



(117) 



^j^Rad<V> I 



■{2S-2P) 
^^|ad<5+y>(25-2P) 
■(25 -2P) 



^j^Rad<SV> I 



1451.7912 MHz 
1089.6513 MHz 
1226.0871 MHz 



(118) 
(119) 
(120) 
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It seems that Eq.(117) is smaller whereas Eq.(118) too large. So as an empirical rule in 



our semiquantitative calculation, we may use Eq.(119) to get 



ltheor.^2S - 2P) = 1089.651 + 0.145 - 0.002 = 1089.794 MHz 



(121) 



which is larger than Eq.(114) by 3%. 



(b) The Lamb shift of He^ atom has been measured as (quoted from j27]): 



L^f+(25 - 2P) = 14041.13(17) MHz 

Similar to the case of hydrogen atom but with Z = 2 and bHe+ 
find 



(122) 
0.0001371, we 



AE^li<^>{2S -2P) 



\TpRad<S+V> 
He+ 

A rpRad<SV> 



{2S - 2P) 
{2S - 2P) 
{2S - 2P) 



1.252680693 x 10^° Hz 

2.023083608 x 10^° Hz 

1.369980830 x 10^° Hz 

1.636521214 x 10^° Hz 



(123) 



As in the case of H atom, we take the < S + V > scheme and add 



AE';^^f-\2S - 2P) = -2.165 kHz (124) 

AE^f+{2S-2P) = 4.514 MHz (125) 

1.2/m), to find the theoretical value: 

L^^IT{2S - 2P) = 13704.220 MHz (126) 



which is smaller than Eq.(122) by 2.41%. 



(c) The following energy-level difference is related to the "hyper Lamb shift (HLS)" [TO] : 



A''^p = Eh{4S)-Eh{2S)-^[Eh{2S)-Eh{1S)]= 4797.338(10) MHz (127) 



Theoretically, now Eq.(102) makes the main contribution: 



25 



AE^^^[(45) - ^(25) + ^(15)] = 3923.95 MHz 



(128) 



(The notation in parenthesis is self-evident). Eq.(103) and Eq.(105) contribute 



^ j^recoil—l 



[(45) - ^{2S) + ^(IS)] = -4.186 MHz 



(129) 



and 



^^Rad<s> ^ 451.229097 MHz 



^^Rad<s+v> ^ 529.288296 MHz 
^^Rad<sv> ^ 675.907131 MHz 
^^Rad<v> ^ 903.266275 MHz (130) 



respectively. Adding a small contribution from Eq.(107) 



AE^^[{AS) - ^{2S) + ^(15)] = 0.1270967854 MHz (131) 



we get 



^Theor.<s> ^ 4371.120197 MHz 

^Theor.<S+V> ^ 4449.I79396 MHz 

^Theor.<sv> ^ 4595.798231 MHz 

^Theore.<v> ^ 3923.95 -4.186 + 903.266275 + 0.1271 = 4823.1574 MHz (132) 



The < V > scheme is only larger than Eq.(127) by 0.54%. All other schemes would be 
too small. So we guess that for S states < V > scheme is better than < S > scheme, 
(d) The following energy-level difference was also measured as [10] : 



A'h" = Eh{AD^/2) - Eh{2S) - ^[Eni'lS) - Eh{IS)] = 6490.144(24) MHz (133) 



Theoretically, Eq.(102) also makes the main contribution: 



AE§^^[{4D5/2) - -(25) + -{IS)] = 5747.92 MHz (134) 
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while 



^ -pprecoil—l 



[{AD^,2) - \{2S) + -^{IS)] = -4.18611 MHz 



(135) 



AE™-'[(4D5/2)] = a^me(5.446170255 x 10 



are all small, we will have 



-6.9283 kHz 



(136) 



^1 j^rad<S> [ 



■[(4/^5/2)- 
■[(4/^5/2)- 
■[(4/^5/2)- 
■[(4/^5/2)- 

Finally, the nucleus size effect gives 



^1 j^rad<S+V> \ 



^1 ^rad<SV> \ 



^(25) + ^(15)] 
\{2S) + \{IS)] 
\{2S) + \{IS)] 
\{2S) + \{IS)] 



302.088631 MHz 



700.843464 MHz 



369.124660 MHz 



500.131264 MHz 



(137) 



A'E^^%AD,/,)~\{2S) + \{IS)] 



In sum, we have 



11.62027752 x 10^( x 

M 4 



0.10894 MHz (138) 



<s> 

H 

i<V> 
H 

><S+V> 
^ H 



A' 
A 



A' 



<sv> 

H 



A'E§^^ + A'E^" 
6444.679 MHz 
6112.961 MHz 
6243.967 MHz 



^/ j^recoil—2 



^/ j^rad<S> 



^1 j^rad<NS> 



6045.925 MHz 



(139) 



which are smaller than the experimental value (133) by 6.8%, 0.7%, 5.8% and 3.8% 
respectively. . 



(e) Experimentally, the combination of Eq.(127) with Eq.(133) yields: 



Then, theoretically, we have 



1692.806 MHz 



(140) 
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A"5^^(4D5/2 - 45) = 1.823886903 x 10^ Hz (141) 

A"^''"*'"^ (4^)5/2 -45) = 1.1008 Hz (142) 

A"^''°*'-'(4D5/2-45) = -6.9283 kHz (143) 

A"2^ (4^5/2 - 45) = -0.0181605862 MHz (144) 



and 



A"™'^^^^ (4^5/2 -45) = -149.1404661 MHz (145) 

A"^°'<^>(4L)5/2 -45) = -202.4228107 Mi/^ (146) 

A"^'^<^+^>(4L)5/2 -45) = -160.1636366 M/7z (147) 

A"™'^<^^>(4D5/2 -45) = -175.7758676 Mifz (148) 



Altogether, we have 



A //theore.<S> 
^ H 

A iitheore.<SV> 
^ H 

//theore.<V> 



(4i^5/2 


-45) 


= 1674.721349 


MHz 


(4/^5/2 


-45) 


= 1663.716339 


MHz 


(4i^5/2 


-45) 


= 1648.104108 


MHz 


(4/^5/2 


-45) 


= 1621.439105 


MHz 



/S.''h {^D^/2-AS) = 1621.439105 MHz (149) 



which are smaller than Eq.(140) by 1.1%, 1.7%, 2.6% and 4.2% respectively, 
(f) It's time to go back to the precision data of 25 — 15 transition in hydrogen atom as 
discussed in section II. Rewrite Eq.([7]) as (see also |13]): 

^E''^^{2S - 15) = 2.46606141318734 x 10^^ Hz (150) 
Theoretically, we have [see Eq.([l7))]: 



AE^^^(25 - 15) = 2.466067984 x 10^^ Hz (151) 



^^recoii-i(^25 - 15) = 22.32598676 MHz (152) 
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^^md<5>(^25- 15) = -5142.081146 MHz 

^£;rad<5+y>(25- 15) = -5765.958928 MHz 

^^md<5V>(^25- 15) = -6835.535314 MHz 

^^rad<y>(25- 15) = -8541.095068 MHz (153) 



AE^^(25 - 15) = 11.62027752 x 10^(- - 1) = -1.016774283 MHz (154) 

8 

radi 



If taking the value of AE™'^(25 - 15), we get 



^^theore.<5>(25- 15) = 2.466062836 x 10^^ Hz 

^^theore.<5+y>(25- 15) = 2.466062239 x 10^^ Hz 

^^tfeeore.<5V>(25- 15) = 2.466061169 x 10^^ Hz 

^^theore.<v>^25- 15) = 2.466059464 X 10^^ Hz (155) 



They are larger than Eq.(150) by 1450 MHz, 826 MHz and smaller than Eq.(150) 
by 244 MHz, 1949 MHz respectively. Or, their discrepancies are +5.9 x 10~^, +3.3 x 
10^^, —1.0 X 10^'', —7.9 X 10^^, respectively. This discrepancy is basically stemming from 
the uncertainty in the calculation of AE^'^{2S — 15). 



(g) Let us turn to the isotope-shift of 25 — 15 transition. Rewrite Eq.([8]) as 

A£;^?^(25 - 15) = 6.70994337 x 10^^ Hz (156) 



Theoretically, rewrite Eqs.(19) and (20) as 



AE^^§{2S - 15) = 6.7101527879 x 10^^ Hz (157) 

and 



AE'j^Zh~\2S - 15) = -11.176 MHz (15^ 
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15) = 


-1.399158 


MHz 






IS) = 


-2.324028 


MHz 






IS) = 


-1.568915 


MHz 






IS) = 


-1.859945 


MHz 


(159) 



AE^^hC^S - IS) = -5.11384949 MHz 



(160) 



Altogether, we find [using <V> scheme in Eq.(159)]: 



A rptheore.<V> 



(25 - IS) = 6.709966701 x 10^' Hz 



(161) 



which is larger than Eq.(156) by 2.333 MHz or only 3.5 x 10 ^. Evidently, even Eq.(157) 
solely deviates from Eq.(156) by 3 x 10^^ only. And as expected, the different schemes for 



AE^_jj{2S — 15*) have little infiuence on the theoretical value, because any one of Eq.(159) 



is much smaller than the nucleus size effect Eq.(160) (r^ = = 2.115/m). 



(h) Finally, the so-called absolute Lamb-shift of IS* state in hydrogen atom was deter- 



mined by Weitz et a/. [TO] from the measured value Eq.(127) or (133). In our notation, using 



Eq.(133), we will write it as follows: 



Lh{1S) = 4{A--^ - AE§^^[{4D,/,) - |(25) + |(15)] - Ai?— '-^[(4^5/2) - |(25) + ^(1^)] 
-AE'h"""'''\4D,/2) + ILh{2S) - LH{iD5/2)} 

(162) 

Here the Lamb shift of 25* state Lh(2S) can be determined from the experimental value 
of Eq.(114) with Lh{'2,Pi/2) being calculated from Eq.(105): 



Lh{2S) = L^^^(25 - 2P) - A£;™-2(25 - 2P1/2) + Lh{2Pi,2) = 1040.901 MHz (163) 



And Lh{4:D^/2) can also be calculated from Eq.(105), so 



Lh{IS) = 8188.478 MHz 



(164) 
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which is in agreement with 8172.874(60) MHz given by [TO] within an accuracy < 0.2%. 
If we use Eq.( |127 ) to derive Lh{1S), we would have to calculate L^i'iS) which is much 
larger than Lh{4:D^I2) and its derivation from Eq.(105) seems not reliable. Similarly, the 
theoretical value of Lh{)-S) turns out to be 



(165) 



with AE^^{IS) = 0.14525347 MHz. However, the value of AE^"'^(1S) strongly depends 



on the scheme we used in Eq.(105), which must be narrowed in a high-loop calculation. The 



theoretical prediction was given in [27j as: 



Itneor.^^g^ = 8172754(14) (32) kHz 
Further discussions can be found in Refs. [SI SU 135] . 



(166) 



IX. SUMMARY AND DISCUSSION 



The remarkable progress of the experimental research on energy-level differences in hy- 
drogenlike atoms has been making this field an ideal theoretical laboratory for physics: 
(a) The inevitable and successful use of reduced Dirac equation (RDE) to hydrogenlike 



atoms, especially to the isotope-shift of 25'— IS transition as reflected by Eqs.(156) through 



(161), is by no means an accidental fortune. It implies that the argument in section III for 



introducing RDE, Eq.(41 ), is correct to a high accuracy. In particular, the basic principle of 



invariance under space-time inversion Eq.(42) (with original mass m) could remain valid even 



for a noninertial frame. This implication has a far-reaching consequence that a generalization 
at the above symmetry to a localized curved space-time may be served as a possible road to 
quantize the general theory of relativity |16j . 

However, there are two realizations of potential V under the space-time inversion. 



Eq.( 35) ("scalar" type) and Eq.( 42) ("vector" type). While Eq.(42) does dominant in an 



atom like H with rrip ^ rrie, the remaining discrepancy of 2.333 MHz between theory and 
experiment [Eq.( 1161 ) versus Eq.(156)] strongly hints that an important and subtle effect had 
been ignored. (To consider the contribution of the deuteron polarizability merely accounts 
for about 20 kHz [6]). We think what neglected must be a tiny excitation of antiparticle 
field in the nucleus due to its interaction with electron in the CMCS. So when we reduce 
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the degrees of freedom of two-body system from two to one, the RDE should be modified 



to take account of the tiny mixture of "scalar" potentials (see the page note after Eq.(42)). 
We don't know how to improve RDE yet. However, an experimental evidence for the above 
conjecture could be the following prediction: The discrepancy between present theory (with 
RDE) and experiment must be smaller for the isotope shift in 25* — 15* transition of atoms 
'^He and ^He than that of atoms H and D. 

Recently, by using Dirac's method, Marsch rigorously solved the hydrogen atom as a two- 
Dirac particle system bound by Coulomb force His solutions are composed of positive 
and negative pairs, corresponding respectively to hydrogen and anti-hydrogen as expected. 
However, surprisingly, in the hydrogen spectrum, besides the normal type-1 solution with 
reduced mass /i, there is another anomalous type— 2 solution with energy levels: E'n = 
Mc^ — 2/ic^ + |yuc^(^)^ + ■ ■ ■ (ra = 1, 2, . . .) and "strange enough, the type— 2 ground state 
(n = 1) does not have lowest energy but the continuum {n = oo)". In our opinion, based 
on what we learnt from the Dirac equation and RDE, these anomalous solutions imply a 
positron moving in the field of proton. So all discrete states with energy E'n are actually 
unbound, they should be and can be ruled out in physics either by the "square integrable 
condition" or the "orthogonality criterion" acting on their rigorous wave functions (for one- 
body Dirac equation, see [35], also p. 28 — 31, 50 of [36]). On the other hand, all continuum 
states {n = oo) with energies lower than Mc^ — 2fic^ correspond to scattering wave functions 
with negative phase shifts , showing the repulsive force between positron and proton, (see 
[37] , section 1.5 in ^B] or section 9.5 of [15j). Marsch's discovery precisely reflects two 
things: (a) the negative energy state of a particle just describes its antiparticle state, (b) 
The Coulomb potential allows a complete set of solutions comprising of two symmetric 
sectors, hydrogen and antihydrogen.In the hydrogen sector, the proton remains unchanged 
regardless of the changing process of electron into positron under the Coulomb interaction. 



The above particle- antiparticle symmetry (including Eq.(42) showing the unequal treat- 
ment between electron and nucleus), together with the parity symmetry, is hidden in the 
Dirac's four-component theory in covariant form so they were overlooked to some extent in 
the past. The advantage or flexibility of two-component noncovariant form of Dirac equation 
or RDE (as discussed in this paper) lies in the fact that the above two symmetries become 
accurate and so easily to be extended (or violated) in an explicit manner. For completeness, 
let us stress again that for antiparticle, one should use the momentum and energy operators 
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being Pc = iV and = versus p = — iV and E = i-^ for particle as required by 
the space-time inversion symmetry. The historical mission of the conception to imagine 
the positron as a "hole" in the sea of negative energy electrons is already over. Since the 
CPT invariance had been further verified [39j, the relation between a particle \a) and its 
antiparticle \d) is well-established as: [§] 



\a) = CPT\a) 

with their wave-functions (in free motion) being respectively: 



(x, t\a) ~ exp[-(p ■ X — Et)] 



(x, t\a) ~ exp[— -(p ■ X — Et)] 

f t 



(167) 



(168) 



(169) 



Note that in Eqs.(168) and (169), they have the same momentum p and positive energy E. 
Either a newly defined space-time inversion (x — )■ — x, t — > —t) or a simple change of i — )■ — i 



will transform Eq.(168) into Eq.(169) (or vice versa). 

(b) Throughout this paper, the electron bound in an atom is just treated like a stationary 
"ball" with nucleus at its center and having a (Bohr) radius (~ l/arrie). However, it is in an 
off-mass-shell state (In some sense, our atom model is just the opposite to J. J. Thomson's 
atom model 100 years ago). In fact, the electron's mass is reduced suddenly from rrie to 
H in the RMCS when it is captured by a nucleus at the far remote orbit with quantum 
number n — > oo and further reduced to + ~ /i— ^^1^ until n decreasing to the lowest 
limit n = 1. The Lamb shift should be viewed as a further modification on the mass of an 
off-mass-shell electron due to radiative correction. 



Notice that the parameter in the vertex function, Eq.(87), means the square of (three- 
dimensional) momentum transfer when a free electron is on its mass-shell and collides with 
some other particle as discussed in Ref . |8bj . By contrast, now exhibits itself as an effective 
potential of 5-function type exerted by the nucleus to the bound (and so off-mass-shell) 



electron as shown by Eq.(92). To bind an electron to a nucleus is a nonperturbative effect. 
Hence we can understand why the discrepancy between (^'^^ (calculated by perturbative 
QED at one-loop order) and C*~^^ (evaluated via nonperturbative Virial theorem) is so 



[§] To our knowledge, the correct definition, Eq.(167), was first given by T. D. Lee and C. S. Wu at Ann. 
Rev. Nucl. Sci. 15, 381(1965). See also G. J. Ni at J. Fudan Univ. (Natural Science) No.3-4, 125(1974). 
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large. Fortunately, they lead to discrepancies in the calculated values of Lamb shift being 
not so large as shown in Section VIII. When (^<^> or (^<^+^> (or (^<^^>^ is substituted 



into the Eq.(105) which is derived from perturbative {L = 1) theory, we should always be 
aware of some theoretical inconsistency in such a semi-empirical treatment. But as a whole, 
we believe that the concept of Lamb shift as an off-mass-shell effect in covariant QED is 
basically correct. 

(c) For a free on-mass-shell electron, its charge square will increase with the increase 



of as shown by Eq.(72) (with /is = me, = —Q"^) and was calculated in detail in [8b], 
coinciding with the experimental data. Note that, however, the Ward identity Zi = Zi had 
been used. An interesting question arises for a bound electron: as its is not a function of 
Q^, will change with the variation of the quantum number n? To answer this question, 
let us put Ward identity aside for a while and write down the renormalized olr = ^ as 

Z^ 

an = ^Z^a — > -fa (170) 
^1 ^1 



Let us work in the CMCS, so Z2 = t~r B is shown in Eq.(53) but with n replaced 



by rrie = m. Similarly, Zi is given by the first part of Eq.(87) with /i — )■ m: 



rv 11 ?T?^ 

^1 ^ 1 + £[y - In ^ - 3^ + 4(1 + In^] (171) 
where the off-mass-shell parameter ^ in CMCS is defined by 



p2 = ^2(^1 - ^) = rn^il -r]-C) = m^{l - v) - ^ C = - ^ C = (1 - C) (172) 



with 



m2(l-r/)=/i^ r/ = l-4, C = 4c 



(173) 



and C is exactly that in Eq.(88). If we ignore the dependence of (1 — B) on C, Eq.( |170 ) 
would give {( <^ 1): 



2a; , , Cm 
Or = a[l + — In 1 + - 

IT 7] 

after renormalization by adjusting the arbitrary constant /ii so that 



(174) 
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aij|f_^o = a (175) 

which connects to the Thomson hmit afjlq^Q = a for a free electron continuously but 
not smoothly. Then for two lowest bound states with n = 1 and n = 2, we would have (in 
<V> scheme): 

aR\n=i ^ a(l. 000433832), aij|„=2 ^ a(1.0001123) (176) 
This would modify the Bohr energy level in hydrogenlike atom A to 

^r^(^) = -^^^A (177) 

and make an extra contribution to the isotope-shift as 

AE^"_^^{2S - IS) ~ 726 MHz (178) 
which is definitely excluded by the experiment. Hence the above consideration from 



Eq.(170) till Eq.(178) is wrong. We learn concretely once again that the Ward identity 
Zi = Zi is valid not only for an electron on its mass-shell, but also for off-mass-shell case. 
Thus we use the same value of a throughout the whole calculation. 

(d) In Ref . |8b] , using our RRM and new renormalization group equation (RGE) for QCD 
derived from it and keeping all masses of 6 quarks [vric = 1.031 GeV, rrib = 4.326 GeV, 
rrit = 175 GeV, iris = 200 MeV, rriu = 8 MeV, rrid = 10 MeV), we calculated the strong 
coupling constant agXQ) foi' i = u,d,s,c,b respectively. Their running curves (starting 
from the common renormalization point as{Mz) = 0.118) follow the trend of experimental 
data (as shown on p. 158 of [39]) quite well but separate at the low Q region. Interesting 
enough, each of them rises to a maximum a™"^ and then suddenly drops to zero at Q = 
corresponding to a threshold energy scale E^^ which could be explained as the excitation 
energy scale for breaking the quark pair. For example, we find E^^ = 1.13 GeV which is 
just the hadronization energy scale of Upsilon T{bb) against its dissociation into two bosons. 
Experimentally, M(T(4s)) - M(T) = 1.12 GeV and T(4s) B+B- or B^B°. similarly, 
= 0.398 GeV while M(V'(3770)) - M(^(3097)) = 673 MeV and ^(3770) ^ D+D" 
or D^D°. it seems that Ef ~ 90 MeV and E^^^^ ~ OA MeV are not so rehable but still 
reasonable. 
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Actually, our calculation on QCD is backed by that on QED. In [8]b, using our RRM 
and improved RGE, keeping all contributions from 9 charged leptons and quarks we were 
able to calculate the running fine-structure constant otR{Q) from the renormalization point 
aR{Q)\Q=Q = a = (137.036)^^ until it coincides with the experimental value of aexpiMz) = 
(128.89)"^. We fitted quark's masses as mentioned above and found no further room left for 
extra charged elementary particles (say, of 4th generation). 

(e) In 1989, we had estimated the upper and lower bounds on Higgs mass by using 
a nonperturbative approach in QFT — the Gaussian effective potential (GEP) method, 
yielding [iO] : 

76 GeV <Mh < 170 GeV (179) 

Like many other authors, we were bothered a lot by divergences. After a deeper study on 
the model by using our new RRM [8a] , we restudied this problem by combination GEP 
with RRM, yielding [231 : 

Mh = 138 GeV (180) 

This is not a upper or lower bound but a prediction based on the input of experimental 
data: 

Mw = 80.359 GeV, Mz = 91.1884 GeV, a'^ = — ^ = 128.89, 



/2 

sin^^vF = — o = 0.2317 
9^ + 9' 



where 6w is the weak mixing (Weinberg) angle. Because of getting rid of all divergences, 
our calculation is well under control at every step. As now the search for Higgs particle 
becomes so urgent experimentally but the theoretical estimation about its mass still remains 



uncertain [39], we think our approach with the prediction (180) deserves to be reconsidered. 

(f) Moreover, our RRM can be used in D + 1 space-time without limitation on the 
space dimension D. A detailed analysis of sinh(sine)-Gordon models with D = 1,2 and 
3 (also using GEPM) is given by Ref.[25]. Another example is again the Lamb shift but 
calculated by QED in noncovariant form and by using RRM similar to that in this paper, 
see Appendix ([26], see also the Appendix 9A in [IS])- The theoretical value (A. 20) seems 



better than (121), showing that for dealing with the Lamb shift, the noncovariant theory 
may be more suitable than the covariant one at least in the lowest order. 

(g) Previously, the theories for Lamb shift or generally for calculating energy levels in 
hydrogenlike atoms are rather complicated as reviewed in refs [H [2] and [27] , some of them 
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have been discussed in the Appendix of this paper. For further clarity, let us try to summarize 
the main obstacles, or challenges in four points: 

(1) Different masses of nuclei must be taken into account; 

(2) Relativistic effects of the electron (not nucleus) are important; 

(3) In calculating radiative corrections, the divergence becomes severer and severer with the 
increase of loop number; 

(4) Since nuclei's properties are different from one atom to another, to treat each atom as a 
two-body system individually would be a daunting task, it couldn't be rigorous eventually 
too. This can be clearly seen from the recent work by Marsch [5^ . 

Facing these challenges and learning from lessons and experiences of previous authors, 
we see that the clue point is to replace the electron mass m by reduced mass fi and work 
in the noninertial frame (RMCS) throughout the entire calculation. As is well known, 
this can be handled in nonrelativistic QM by a mathematical trick but is impossible in 
relativistic case. So what we need is a new understanding on the essence of special relativity 
— the invariance (of theory) under the (newly defined) space-time inversion in one inertial 
frame. Then we are able to claim the same invariance in the RMCS with /i replacing m for 
establishing the RDE, ignoring a small centripetal acceleration of the nucleus in CMCS (see 



page note after Eq.(42)). The approximation in RDE is some price paid for the much bigger 
gain — improving the original Dirac equation (unable to treat different nuclei) and avoiding 
the confusion in QED calculation because of the entanglement of two frames: CMCS with 
RMCS ( i.e., the radiative corrections are entangled with the recoil effect as we can see from 
previous literatures). In some senses, we jump over obstacles (1) ,(2) and (4)at the least 
labor cost (by constructing RDE). In the meantime, we hope RDE would help to ease the 
difficulty in point (3). And it's a great pleasure to see that the essential correctness of our 
understanding has been validated by Marsch's work as well as puzzles raised in his paper 
[31]. Please see also Ref [IT] . 

As to challenge (3), only after we puzzled over the "divergence" for decades, could we 
suddenly realize that we misread its implication as a " large number" . Rather, it means the 
"uncertainty". Let us look at the calculation in section IV again. Previously, many authors 



treated the divergent integral I in Eq.(50) by different tricks of regularization , arriving at 
Eq.(56). Because both A and B are divergent, it was thought that the original mass (// 
here) does receive some radiative corrections (via the self-energy diagrams in Fig 2(a) and 
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(b)) and becomes a "renormalized" mass {fiR here). The latter should be the observed mass 
in experiment or physical mass (of electron). So the original mass was called as the "bare 
mass", which was written into the Lagrangian density as an input parameter of QFT. Then 
in constructing Feynman diagrams of certain perturbative calculation, one needs to further 
introduce some (divergent) "counter terms" for cancelling the divergence stemming from 
the bare mass. Based on that understanding, the renormalization factor for wavefunction. 



Z2 in Eq.(55), would be a divergent quantity too (in sharp contrast to here Eq.(59) being 



a fixed and finite number). Previously, In Eq.(72), while the on the left handed side is 
the observed charge which should be finite, the e on the right handed side was regarded as 
a "bare charge" which, together with the Z3, was a divergent quantity, (see Fig. 7.8 in |23] . 
By contrast, here both Z3 and e are finite. Actually, here e is defined as the physical charge 
observed at the Thomson limit in experiment). 

In our opinion, the reason why we encountered so many superfluous troubles in the 
past is because we overlooked what Bethe said in 1947(22] ■ Please read his words quoted 



after Eq. (A. 2) in the Appendix. Let us explain via our Eq.(46). The (reduced) mass 
(/i) already contains some contributions from self-energy diagrams like Fig. 2(a) and (b). 



When we evaluate the (divergent) integral, Eq.(50), trying to find the radiative corrections 
on the electron, the latter is bound to confuse with that already contained in the mass. 
In other words, the dividing line between them is blurred inevitably. The emergence of 
explicit divergence is essentially a warning that the effect you want to evaluate has been 
entangled with the mass itself, rendering both of them uncertain. Hence the aim of so-called 
renormalization is nothing but to redraw the dividing line between them such that the values 
of mass (reconfirmed by the experiment) and the new effect {e.g., the mass increment when 



the electron is moving off-mass-shell, Eq.(61)) can be clearly separated. In short, what we 
have been learning in the past decade is: At the level of QM, in the Hamiltonian like Eq. 
(A.l), the parameters m and e can be regarded as well-defined. But they are not so at the 
level of QFT. Once the calculation is made beyond the tree level, i.e., with loop number 
L > 1, the divergence occurs and the meaning of parameters becomes obscure immediately. 

We need to reconfirm all parameters contained in the Lagrangian density before they can 
be linked with experiments. In this sense, a model of QFT is at most an "effective field 
theory" . According to the above point of view, we believe that our RRM just provides a 
natural way to carry out these processes of reconfirmation [Saj, getting rid of divergences 
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and ambiguities. Please see also Ref |i2] . 

(h) Last, but not least, during the learning and teaching of graduate courses on QFT 
for decades, we have been sharing the joy and puzzle with our students all the time. We 
hope that the presentation of this paper could be useful as a teaching reference to render 
the QFT course more understandable, interesting and attractive. 
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Appendix: Comparison Between Noncovariant and Covariant Theories for Lamb 

Shift 

1. To our knowledge, the precision theory for Lamb shift was based on a combination 
of noncovariant (nonrelativistic or old-fashioned) QED with covariant (or relativistic) QED 
as discussed in Ref.pTj- As explained clearly by Sakurai in Ref.pl], in perturbative QFT 
of noncovariant form, all virtual particles are "on-mass-shell". Here we wish to emphasize 
that a rigorous reconfirmation procedure of mass parameter was often overlooked in previous 
literatures. The theory for hydrogenlike atom begins with a Hamiltonian: 

Ho = —p' + - (Al) 

2m ZTTiN r 

(p = —iV, see Eq.(34) in i27j|). As Bethe |29j first pointed out that the effect of electron's 
interaction with the vector potential A of radiation field (see |21j.|15]) 

= — A ■ p {A.2) 

mc 

should properly be regarded as already included in the observed mass rriobs of the electron, 
which is denoted by m in (Al). However, once a concrete calculation is made with (A2) 
being taken into account, the divergence emerges immediately. What does it mean? Math- 
ematicians teach us that there are three implications for a divergence: 
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(a) It is a dimensionless number; (b)It is a large number; (c)It is uncertain. While we physi- 
cists often emphasized the point (b), we didn't pay enough attention to the points (a) and 
(c). We often talked about a quadratically (or linearly) divergent integral without noticing 
that it has a dimension (say, mass dimension) and thus meaningless in mathematics unless 
a mass parameter (say, m) in the integral is already fixed as a mass " unit" so that the inte- 
gral can be divided by (or m) to become dimensionless. Alternatively, a logarithmically 



divergent integral is dimensionless and thus unaffected by the choice of unit [like Eq.(50), 
see also Eq.(A6) below], it just implies an uncertainty waiting to be fixed. The implication 
of uncertainty of a divergence will never vanish even after we introduced a cutoff by hand 
to curb it. For example, in a pioneering paper to explain the Lamb shift, Welton ([30], see 
section 9.6B in Ref.|15j) encountered an integral / = J^'""'' ^ with to being the (angular) 
frequency of virtual photon (vacuum fluctuation). He simply set Umin ~ mZa = Z/a, (a 
is Bohr radius) and Umax ^ m so that / ^ \n{l/Za) = 4.92 (for Z = 1) which leads to an 
estimation of Lamb shift L^'^"'"- (25*1/2 — 2P1/2) — 668 MHz. If instead of Bohr radius, the 
lower cutoff is provided by the electron binding energy, one should get I ~ In(Za)^^ and 
jjheor. i33g MHz. (scc Eq.(30) in [27]). The above arbitrariness just reflects what essen- 
tial in a divergent integral is not its large magnitude {\a{Za)~^ is merely of the order of 10) 
but its uncertainty. So what important in handling the integral is not to curb (or to hide) 
its divergence but let the divergence exhibits itself as some arbitrary constants explicitly (as 
shown in section IV- VI). We will show later how to do this way for noncovariant QED. 
2. While Eqs.(Al) and (A2) only describe a spinless particle, the electron has spin which 



endows it with the relativistic nature as shown by Eqs.(41 )-(45). For two-particle system, 
based on Bethe-Salpeter equation, an effective Dirac equation (EDE) was derived as shown 
by Eq.(23) in [27|. When the electromagnetic field interaction is taken into account, the Breit 
potential Vbt was derived as shown by Eq.(35) in [27J. Then the total Breit Hamiltonian 
reads (Eq.(36) in [27]): 

HBr = Ho + Vsr (A.3) 

However, the electron mass m' (in our notation here) appeared in EDE or Vbt should be 
that in the Dirac equation, also that in the definition of reduced mass = , eventually 

m' could be identified with the observed mass rriobs, which is not equal to the m in Eq.(Al). 
This is because besides (A2) there is an extra interaction due to electron spin with the 
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radiation field: 

{g = 2 X 1.0011596522 is the gyromagnetic ratio of electron, see Eq.(9A.15) in [U]). The 
difference between m and m' will be calculated in (A16) below. It turns out to be of the 
order of am and cannot be ignored at the level of QED, especially for the explanation of 
Lamb shift. We guess this must be one of the reasons why all calculations based on Eq.(A3) 
became so complicated. 

3. In noncovariant theory, the leading contribution to the Lamb shift comes from the one- 
photon electron self-energy. The nomenclature here is different from that in the covariant 
theory. Roughly speaking, so-called electron self-energy often corresponds to the vertex 
function in covariant theory (Fig.2(d) in this paper) or to Figs. 8 and 11 in Ref. [27] and its 
evaluations have extended over 50 years [31] . More precisely, it is identified with the radiative 
insertions in the electron line and the Dirac form factor contribution. Further contributions 
from the Pauli form factor and the vacuum polarization [27] will add to a theoretical value 
of classic Lamb shift being 1050.559 MHz. If taking more high-order corrections into ac- 
count, the theoretical value coincides with the experimental value 1057.845 MHz rather 
accurately (see Table 20 in [27]). However, the above calculation looks quite complicated 
due to two reasons: (a) The difficulty of dealing with two masses in two coordinate systems, 
the electron mass m and the reduced mass fi; (b) The introduction of an auxiliary parame- 
ter a [m{Za) ^ a ^ m{ZaY] to separate the radiative photon integration region into two 
parts. In the low momentum region, the Bethe Logarithm [32j in noncovariant form makes 
the main contribution. In the high momentum region, the evaluation is resorting to some 
relativistic covariant form [22] . Then two expressions are matched together to get the correct 
result. It seems to us that the matching trick used is doubtful because both ultraviolet and 
infrared divergences were ambiguously handled by some cutoff which missed the main point 
of renormalization — to reconfirm the mass parameter in the presence of radiative corrections 
as shown in section IV (covariant form) or below. 

4. A simple calculation for Lamb shift in noncovariant form was proposed in Ref. |26| (see 
also Appendix 9A of Ref. [15] ) . Consider the self-energy diagram of an electron with reduced 
mass and (three-dimensional) momentum p in the RMCS of a hydrogenlike atom. Similar 
to Fig.2(a), but also different in the virtual state, now a photon has energy ouk = k = |k| 
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while the electron has momentum q = p — k and energy Eg = ^(f- The electron in plane- 
wave state |p > has two interactions with the radiative field at each vertex as shown by 
(A2) (m — )■ /i) and (A4), acquiring an increase in energy respectively (see FIG. 3): 

A^O-) = ^ l<^l^.^.|P>l ^ (, = 1,2) (A5) 

Here Ei = Eg + is the energy of the intermediate virtual state \i >. Simple evaluation 
leads to 

where 77 = cos 6* with 6 being the angle between k and p, = 2(yU — prf). Like Eq.(50), we 
take partial derivative of the divergent integral / with respect to ^ (then the integration of 
k) and integrate back to / again, yielding: 

AE« = feS'V + + • • • (A.7) 



n^ a 4 4 4 

Note that the term Ij^^p^ will combine with the kinetic energy of a ("spinless") electron, 
they are indistinguishable. The appearance of an arbitrary constant Ci precisely refiects the 
fact that we cannot find the reduced mass via the valuation of AEp^^ in perturbation theory. 
So we must choose b^i^ = to reconfirm the value of /i (which is still not the final observed 
mass, see below). Similar evaluation on if^-^] (of the real electron with spin 1/2) which would 
induce the spin fiip process between states |p, ±^ > and |q, ±| >, leads to 

J= / (A.IO) 







Being a quadratically divergent integral, J needs partial derivative of third order with respect 
to yielding: 

AEf = 6(') + 6SV + 6(V + --- (All) 
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if = l"Ji\iil^2 + - 4C, - - ^1 (^.12) 

4 TT IjL jJ,"^ 



io\ a'^ a A 4 4 

= T — + 2 + - In/i - {A.13) 

^ 4 7r/i3^ 15^ ^ ^ 

Let's manage to fix three arbitrary constants C2, C3 and C4. First, the term bf'^'p'^ should be 
combined with term. Since fi is aheady fixed, further modification on /i due to electron 
spin should be finite and fixed. So the only possible choice of C2 is to cancel In fi which is 
ambiguous in dimension: C2 = In fi, yielding 

f 2) 

Then the dimensional constants C3 and C4 must be chosen such that % =0, implying that 
the starting point of this theory is the nonrelativistic Hamiltonian Hq in Eq.(Al) without 
rest energy term while both masses of the nucleus and the electron (with spin) are fixed by 
experiments. Hence now /i acquires a modification via b^^^p"^ term and becomes an observable 
one: 

/i — > fiobs = {A.16) 



However, we have to consider the relativistic energy of electron shown in Eq.(30), where 
the term (— g^P^) goes beyond Eq.(Al). Yet the modification of fi shown as (A16) does 
induce a corresponding change — h)p'^y which should be regarded as an invisible 

^ obs ^ 

"background" and subtracted from the term induced by radiative corrections. (The 
relativistic correction is brought in via the RDE as discussed in section VIII). As a whole, 



the combination of contributions from if^-J and HlJ^ leads to 



, (1) , , (2) _ , (2) 



61 = bY> + bY> = b\'> (A.17) 
and a "renormalized" 62: 

6« = 6« + 6f ) + ^(3/3 + 3/32 + (3^) - ^(1.99808) (A.18) 
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Here we only keep the lowest approximation at the last step. Hence the electron self-energy- 
diagram contributes a radiative correction to the energy level of the stationary state \Z,n,l > 
in a hydrogenlike atom: 

AE^^'{Z,n,l) = {Z,n,l\b^p'\Z,n,l) = -3]^^^/^L (^-19) 

This form, together with contributions from the vacuum polarization and nuclear size effect, 
gives a theoretical value for classic Lamb shift: 

^theor.^25i/2 - 2P1/2) ~ 1056.52 MHz {A.20) 

which is smaller than the experimental value by 0.13%. Despite its approximation involved, 
the above method clearly shows that so-called renormalization is nothing but a reconfirma- 
tion process of mass. We must reconfirm the mass before it could be modified via radiative 
corrections. Either "skipping over the first step" or "combining two steps into one" is not 
allowed. 

5. In noncovariant theory, the (three-dimensional) momentum p is combined with the re- 
duced mass /i to form a kinetic energy term j^p"^ on the mass shell. Once the energy is 
modified whereas p is conserved at the vertex, /i is bound to be modified. On the other 
hand, in covariant theory, the electron energy turns to a component of four-dimensional 
momentum p and the latter is conserved at the vertex. So the (reduced) mass /i cannot be 
modified on the mass shell (p^ = yU^). Therefore, the renormalization as some reconfirmation 
has different meaning in covariant theory versus that in noncovariant theory. We guess this 
is why the matching procedure of these two formalisms into one theory for Lamb shift proves 
so difficult. 

6. Every theory in physics is not only a discovery of natural law, but also an invention 
of human being [33]. Hence the comparison among various theories, in many cases, is not 
about a problem of being right or wrong. Rather, it's about a choice of simplicity, harmony 
(self-consistency) and beauty. Only time can tell. 
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FIG. 1: A hydrogenlike atom in quantum mechanical description. The nucleus with mass m2 
occupies a small sphere with radius ttv (greatly exaggerated in the diagram) while the electron 
with mass mi spreads over a larger sphere with radius Re (i.e. atomic radius). Their common 
center is the atom's center of mass (CM). The wavefunction ^/^(r)e~*'^* with r = ri — r2 shows the 
electron's amplitude under a "fictitious measurement" [15], during which the electron and nucleus 
shrink into two "fictitious point particles " located at ri and r2 simultaneously. The Coulomb 
potential V{r) = —^f- between them is a static one. The probability to find the electron at r is 
|'0(r)p while that to find its momentum being p is |i;^(p)p with (j){p) being the Fourier transform 
of ip{r). 



47 




FIG. 2: Four Feynman diagrams at one- loop level (in covariant form), (a) and (b) are self-energy 
diagrams of the electron, (c) is vacuum polarization, (d) is vertex function. Solid lines and wavy 
lines refer to electron and photon respectively, while X denotes the nucleus. Here p, q and k are 
four-dimensional momenta. 




FIG. 3: The electron self-energy (radiative correction) diagram at one-loop level of perturbative 
QCD in noncovariant form. (A. 2) or H^^^ (-^-.4) is inserted into two vertices. Here p, q and 

k are three dimensional momenta. 



